\(\int \frac {x^3 \cosh (c+d x)}{(a+b x^3)^2} \, dx\) [102]

   Optimal result
   Rubi [A] (verified)
   Mathematica [C] (verified)
   Maple [C] (warning: unable to verify)
   Fricas [B] (verification not implemented)
   Sympy [F(-1)]
   Maxima [F]
   Giac [F]
   Mupad [F(-1)]

Optimal result

Integrand size = 19, antiderivative size = 718 \[ \int \frac {x^3 \cosh (c+d x)}{\left (a+b x^3\right )^2} \, dx=-\frac {x \cosh (c+d x)}{3 b \left (a+b x^3\right )}-\frac {\sqrt [3]{-1} \cosh \left (c+\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Chi}\left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right )}{9 a^{2/3} b^{4/3}}+\frac {(-1)^{2/3} \cosh \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Chi}\left (-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right )}{9 a^{2/3} b^{4/3}}+\frac {\cosh \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Chi}\left (\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}+d x\right )}{9 a^{2/3} b^{4/3}}-\frac {d \text {Chi}\left (\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}+d x\right ) \sinh \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{9 \sqrt [3]{a} b^{5/3}}-\frac {(-1)^{2/3} d \text {Chi}\left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right ) \sinh \left (c+\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{9 \sqrt [3]{a} b^{5/3}}+\frac {\sqrt [3]{-1} d \text {Chi}\left (-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right ) \sinh \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{9 \sqrt [3]{a} b^{5/3}}+\frac {(-1)^{2/3} d \cosh \left (c+\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Shi}\left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right )}{9 \sqrt [3]{a} b^{5/3}}+\frac {\sqrt [3]{-1} \sinh \left (c+\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Shi}\left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right )}{9 a^{2/3} b^{4/3}}-\frac {d \cosh \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Shi}\left (\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}+d x\right )}{9 \sqrt [3]{a} b^{5/3}}+\frac {\sinh \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Shi}\left (\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}+d x\right )}{9 a^{2/3} b^{4/3}}+\frac {\sqrt [3]{-1} d \cosh \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Shi}\left (\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}+d x\right )}{9 \sqrt [3]{a} b^{5/3}}+\frac {(-1)^{2/3} \sinh \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Shi}\left (\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}+d x\right )}{9 a^{2/3} b^{4/3}} \]

[Out]

1/9*Chi(a^(1/3)*d/b^(1/3)+d*x)*cosh(c-a^(1/3)*d/b^(1/3))/a^(2/3)/b^(4/3)-1/9*(-1)^(1/3)*Chi((-1)^(1/3)*a^(1/3)
*d/b^(1/3)-d*x)*cosh(c+(-1)^(1/3)*a^(1/3)*d/b^(1/3))/a^(2/3)/b^(4/3)+1/9*(-1)^(2/3)*Chi(-(-1)^(2/3)*a^(1/3)*d/
b^(1/3)-d*x)*cosh(c-(-1)^(2/3)*a^(1/3)*d/b^(1/3))/a^(2/3)/b^(4/3)-1/3*x*cosh(d*x+c)/b/(b*x^3+a)-1/9*(-1)^(2/3)
*d*cosh(c+(-1)^(1/3)*a^(1/3)*d/b^(1/3))*Shi(-(-1)^(1/3)*a^(1/3)*d/b^(1/3)+d*x)/a^(1/3)/b^(5/3)-1/9*d*cosh(c-a^
(1/3)*d/b^(1/3))*Shi(a^(1/3)*d/b^(1/3)+d*x)/a^(1/3)/b^(5/3)+1/9*(-1)^(1/3)*d*cosh(c-(-1)^(2/3)*a^(1/3)*d/b^(1/
3))*Shi((-1)^(2/3)*a^(1/3)*d/b^(1/3)+d*x)/a^(1/3)/b^(5/3)-1/9*d*Chi(a^(1/3)*d/b^(1/3)+d*x)*sinh(c-a^(1/3)*d/b^
(1/3))/a^(1/3)/b^(5/3)+1/9*Shi(a^(1/3)*d/b^(1/3)+d*x)*sinh(c-a^(1/3)*d/b^(1/3))/a^(2/3)/b^(4/3)-1/9*(-1)^(2/3)
*d*Chi((-1)^(1/3)*a^(1/3)*d/b^(1/3)-d*x)*sinh(c+(-1)^(1/3)*a^(1/3)*d/b^(1/3))/a^(1/3)/b^(5/3)-1/9*(-1)^(1/3)*S
hi(-(-1)^(1/3)*a^(1/3)*d/b^(1/3)+d*x)*sinh(c+(-1)^(1/3)*a^(1/3)*d/b^(1/3))/a^(2/3)/b^(4/3)+1/9*(-1)^(1/3)*d*Ch
i(-(-1)^(2/3)*a^(1/3)*d/b^(1/3)-d*x)*sinh(c-(-1)^(2/3)*a^(1/3)*d/b^(1/3))/a^(1/3)/b^(5/3)+1/9*(-1)^(2/3)*Shi((
-1)^(2/3)*a^(1/3)*d/b^(1/3)+d*x)*sinh(c-(-1)^(2/3)*a^(1/3)*d/b^(1/3))/a^(2/3)/b^(4/3)

Rubi [A] (verified)

Time = 0.77 (sec) , antiderivative size = 718, normalized size of antiderivative = 1.00, number of steps used = 23, number of rules used = 6, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.316, Rules used = {5399, 5389, 3384, 3379, 3382, 5400} \[ \int \frac {x^3 \cosh (c+d x)}{\left (a+b x^3\right )^2} \, dx=-\frac {\sqrt [3]{-1} \cosh \left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}+c\right ) \text {Chi}\left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right )}{9 a^{2/3} b^{4/3}}+\frac {(-1)^{2/3} \cosh \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Chi}\left (-x d-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{9 a^{2/3} b^{4/3}}+\frac {\cosh \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Chi}\left (x d+\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{9 a^{2/3} b^{4/3}}+\frac {\sqrt [3]{-1} \sinh \left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}+c\right ) \text {Shi}\left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right )}{9 a^{2/3} b^{4/3}}+\frac {\sinh \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Shi}\left (x d+\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{9 a^{2/3} b^{4/3}}+\frac {(-1)^{2/3} \sinh \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Shi}\left (x d+\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{9 a^{2/3} b^{4/3}}-\frac {d \sinh \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Chi}\left (x d+\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{9 \sqrt [3]{a} b^{5/3}}-\frac {(-1)^{2/3} d \sinh \left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}+c\right ) \text {Chi}\left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right )}{9 \sqrt [3]{a} b^{5/3}}+\frac {\sqrt [3]{-1} d \sinh \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Chi}\left (-x d-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{9 \sqrt [3]{a} b^{5/3}}+\frac {(-1)^{2/3} d \cosh \left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}+c\right ) \text {Shi}\left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right )}{9 \sqrt [3]{a} b^{5/3}}-\frac {d \cosh \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Shi}\left (x d+\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{9 \sqrt [3]{a} b^{5/3}}+\frac {\sqrt [3]{-1} d \cosh \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Shi}\left (x d+\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{9 \sqrt [3]{a} b^{5/3}}-\frac {x \cosh (c+d x)}{3 b \left (a+b x^3\right )} \]

[In]

Int[(x^3*Cosh[c + d*x])/(a + b*x^3)^2,x]

[Out]

-1/3*(x*Cosh[c + d*x])/(b*(a + b*x^3)) - ((-1)^(1/3)*Cosh[c + ((-1)^(1/3)*a^(1/3)*d)/b^(1/3)]*CoshIntegral[((-
1)^(1/3)*a^(1/3)*d)/b^(1/3) - d*x])/(9*a^(2/3)*b^(4/3)) + ((-1)^(2/3)*Cosh[c - ((-1)^(2/3)*a^(1/3)*d)/b^(1/3)]
*CoshIntegral[-(((-1)^(2/3)*a^(1/3)*d)/b^(1/3)) - d*x])/(9*a^(2/3)*b^(4/3)) + (Cosh[c - (a^(1/3)*d)/b^(1/3)]*C
oshIntegral[(a^(1/3)*d)/b^(1/3) + d*x])/(9*a^(2/3)*b^(4/3)) - (d*CoshIntegral[(a^(1/3)*d)/b^(1/3) + d*x]*Sinh[
c - (a^(1/3)*d)/b^(1/3)])/(9*a^(1/3)*b^(5/3)) - ((-1)^(2/3)*d*CoshIntegral[((-1)^(1/3)*a^(1/3)*d)/b^(1/3) - d*
x]*Sinh[c + ((-1)^(1/3)*a^(1/3)*d)/b^(1/3)])/(9*a^(1/3)*b^(5/3)) + ((-1)^(1/3)*d*CoshIntegral[-(((-1)^(2/3)*a^
(1/3)*d)/b^(1/3)) - d*x]*Sinh[c - ((-1)^(2/3)*a^(1/3)*d)/b^(1/3)])/(9*a^(1/3)*b^(5/3)) + ((-1)^(2/3)*d*Cosh[c
+ ((-1)^(1/3)*a^(1/3)*d)/b^(1/3)]*SinhIntegral[((-1)^(1/3)*a^(1/3)*d)/b^(1/3) - d*x])/(9*a^(1/3)*b^(5/3)) + ((
-1)^(1/3)*Sinh[c + ((-1)^(1/3)*a^(1/3)*d)/b^(1/3)]*SinhIntegral[((-1)^(1/3)*a^(1/3)*d)/b^(1/3) - d*x])/(9*a^(2
/3)*b^(4/3)) - (d*Cosh[c - (a^(1/3)*d)/b^(1/3)]*SinhIntegral[(a^(1/3)*d)/b^(1/3) + d*x])/(9*a^(1/3)*b^(5/3)) +
 (Sinh[c - (a^(1/3)*d)/b^(1/3)]*SinhIntegral[(a^(1/3)*d)/b^(1/3) + d*x])/(9*a^(2/3)*b^(4/3)) + ((-1)^(1/3)*d*C
osh[c - ((-1)^(2/3)*a^(1/3)*d)/b^(1/3)]*SinhIntegral[((-1)^(2/3)*a^(1/3)*d)/b^(1/3) + d*x])/(9*a^(1/3)*b^(5/3)
) + ((-1)^(2/3)*Sinh[c - ((-1)^(2/3)*a^(1/3)*d)/b^(1/3)]*SinhIntegral[((-1)^(2/3)*a^(1/3)*d)/b^(1/3) + d*x])/(
9*a^(2/3)*b^(4/3))

Rule 3379

Int[sin[(e_.) + (Complex[0, fz_])*(f_.)*(x_)]/((c_.) + (d_.)*(x_)), x_Symbol] :> Simp[I*(SinhIntegral[c*f*(fz/
d) + f*fz*x]/d), x] /; FreeQ[{c, d, e, f, fz}, x] && EqQ[d*e - c*f*fz*I, 0]

Rule 3382

Int[sin[(e_.) + (Complex[0, fz_])*(f_.)*(x_)]/((c_.) + (d_.)*(x_)), x_Symbol] :> Simp[CoshIntegral[c*f*(fz/d)
+ f*fz*x]/d, x] /; FreeQ[{c, d, e, f, fz}, x] && EqQ[d*(e - Pi/2) - c*f*fz*I, 0]

Rule 3384

Int[sin[(e_.) + (f_.)*(x_)]/((c_.) + (d_.)*(x_)), x_Symbol] :> Dist[Cos[(d*e - c*f)/d], Int[Sin[c*(f/d) + f*x]
/(c + d*x), x], x] + Dist[Sin[(d*e - c*f)/d], Int[Cos[c*(f/d) + f*x]/(c + d*x), x], x] /; FreeQ[{c, d, e, f},
x] && NeQ[d*e - c*f, 0]

Rule 5389

Int[Cosh[(c_.) + (d_.)*(x_)]*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Int[ExpandIntegrand[Cosh[c + d*x], (a
 + b*x^n)^p, x], x] /; FreeQ[{a, b, c, d}, x] && ILtQ[p, 0] && IGtQ[n, 0] && (EqQ[n, 2] || EqQ[p, -1])

Rule 5399

Int[Cosh[(c_.) + (d_.)*(x_)]*(x_)^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Simp[x^(m - n + 1)*(a + b*
x^n)^(p + 1)*(Cosh[c + d*x]/(b*n*(p + 1))), x] + (-Dist[(m - n + 1)/(b*n*(p + 1)), Int[x^(m - n)*(a + b*x^n)^(
p + 1)*Cosh[c + d*x], x], x] - Dist[d/(b*n*(p + 1)), Int[x^(m - n + 1)*(a + b*x^n)^(p + 1)*Sinh[c + d*x], x],
x]) /; FreeQ[{a, b, c, d}, x] && ILtQ[p, -1] && IGtQ[n, 0] && RationalQ[m] && (GtQ[m - n + 1, 0] || GtQ[n, 2])

Rule 5400

Int[(x_)^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_)*Sinh[(c_.) + (d_.)*(x_)], x_Symbol] :> Int[ExpandIntegrand[Sinh[c
 + d*x], x^m*(a + b*x^n)^p, x], x] /; FreeQ[{a, b, c, d}, x] && ILtQ[p, 0] && IntegerQ[m] && IGtQ[n, 0] && (Eq
Q[n, 2] || EqQ[p, -1])

Rubi steps \begin{align*} \text {integral}& = -\frac {x \cosh (c+d x)}{3 b \left (a+b x^3\right )}+\frac {\int \frac {\cosh (c+d x)}{a+b x^3} \, dx}{3 b}+\frac {d \int \frac {x \sinh (c+d x)}{a+b x^3} \, dx}{3 b} \\ & = -\frac {x \cosh (c+d x)}{3 b \left (a+b x^3\right )}+\frac {\int \left (-\frac {\cosh (c+d x)}{3 a^{2/3} \left (-\sqrt [3]{a}-\sqrt [3]{b} x\right )}-\frac {\cosh (c+d x)}{3 a^{2/3} \left (-\sqrt [3]{a}+\sqrt [3]{-1} \sqrt [3]{b} x\right )}-\frac {\cosh (c+d x)}{3 a^{2/3} \left (-\sqrt [3]{a}-(-1)^{2/3} \sqrt [3]{b} x\right )}\right ) \, dx}{3 b}+\frac {d \int \left (-\frac {\sinh (c+d x)}{3 \sqrt [3]{a} \sqrt [3]{b} \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}-\frac {(-1)^{2/3} \sinh (c+d x)}{3 \sqrt [3]{a} \sqrt [3]{b} \left (\sqrt [3]{a}-\sqrt [3]{-1} \sqrt [3]{b} x\right )}+\frac {\sqrt [3]{-1} \sinh (c+d x)}{3 \sqrt [3]{a} \sqrt [3]{b} \left (\sqrt [3]{a}+(-1)^{2/3} \sqrt [3]{b} x\right )}\right ) \, dx}{3 b} \\ & = -\frac {x \cosh (c+d x)}{3 b \left (a+b x^3\right )}-\frac {\int \frac {\cosh (c+d x)}{-\sqrt [3]{a}-\sqrt [3]{b} x} \, dx}{9 a^{2/3} b}-\frac {\int \frac {\cosh (c+d x)}{-\sqrt [3]{a}+\sqrt [3]{-1} \sqrt [3]{b} x} \, dx}{9 a^{2/3} b}-\frac {\int \frac {\cosh (c+d x)}{-\sqrt [3]{a}-(-1)^{2/3} \sqrt [3]{b} x} \, dx}{9 a^{2/3} b}-\frac {d \int \frac {\sinh (c+d x)}{\sqrt [3]{a}+\sqrt [3]{b} x} \, dx}{9 \sqrt [3]{a} b^{4/3}}+\frac {\left (\sqrt [3]{-1} d\right ) \int \frac {\sinh (c+d x)}{\sqrt [3]{a}+(-1)^{2/3} \sqrt [3]{b} x} \, dx}{9 \sqrt [3]{a} b^{4/3}}-\frac {\left ((-1)^{2/3} d\right ) \int \frac {\sinh (c+d x)}{\sqrt [3]{a}-\sqrt [3]{-1} \sqrt [3]{b} x} \, dx}{9 \sqrt [3]{a} b^{4/3}} \\ & = -\frac {x \cosh (c+d x)}{3 b \left (a+b x^3\right )}-\frac {\cosh \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \int \frac {\cosh \left (\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}+d x\right )}{-\sqrt [3]{a}-\sqrt [3]{b} x} \, dx}{9 a^{2/3} b}-\frac {\left (d \cosh \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right )\right ) \int \frac {\sinh \left (\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}+d x\right )}{\sqrt [3]{a}+\sqrt [3]{b} x} \, dx}{9 \sqrt [3]{a} b^{4/3}}-\frac {\cosh \left (c+\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \int \frac {\cos \left (\frac {(-1)^{5/6} \sqrt [3]{a} d}{\sqrt [3]{b}}-i d x\right )}{-\sqrt [3]{a}-(-1)^{2/3} \sqrt [3]{b} x} \, dx}{9 a^{2/3} b}+\frac {\left ((-1)^{5/6} d \cosh \left (c+\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}\right )\right ) \int \frac {\sin \left (\frac {(-1)^{5/6} \sqrt [3]{a} d}{\sqrt [3]{b}}-i d x\right )}{\sqrt [3]{a}+(-1)^{2/3} \sqrt [3]{b} x} \, dx}{9 \sqrt [3]{a} b^{4/3}}-\frac {\cosh \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \int \frac {\cos \left (\frac {\sqrt [6]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-i d x\right )}{-\sqrt [3]{a}+\sqrt [3]{-1} \sqrt [3]{b} x} \, dx}{9 a^{2/3} b}+\frac {\left (\sqrt [6]{-1} d \cosh \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right )\right ) \int \frac {\sin \left (\frac {\sqrt [6]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-i d x\right )}{\sqrt [3]{a}-\sqrt [3]{-1} \sqrt [3]{b} x} \, dx}{9 \sqrt [3]{a} b^{4/3}}-\frac {\sinh \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \int \frac {\sinh \left (\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}+d x\right )}{-\sqrt [3]{a}-\sqrt [3]{b} x} \, dx}{9 a^{2/3} b}-\frac {\left (d \sinh \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right )\right ) \int \frac {\cosh \left (\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}+d x\right )}{\sqrt [3]{a}+\sqrt [3]{b} x} \, dx}{9 \sqrt [3]{a} b^{4/3}}-\frac {\left (i \sinh \left (c+\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}\right )\right ) \int \frac {\sin \left (\frac {(-1)^{5/6} \sqrt [3]{a} d}{\sqrt [3]{b}}-i d x\right )}{-\sqrt [3]{a}-(-1)^{2/3} \sqrt [3]{b} x} \, dx}{9 a^{2/3} b}+\frac {\left (\sqrt [3]{-1} d \sinh \left (c+\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}\right )\right ) \int \frac {\cos \left (\frac {(-1)^{5/6} \sqrt [3]{a} d}{\sqrt [3]{b}}-i d x\right )}{\sqrt [3]{a}+(-1)^{2/3} \sqrt [3]{b} x} \, dx}{9 \sqrt [3]{a} b^{4/3}}-\frac {\left (i \sinh \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right )\right ) \int \frac {\sin \left (\frac {\sqrt [6]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-i d x\right )}{-\sqrt [3]{a}+\sqrt [3]{-1} \sqrt [3]{b} x} \, dx}{9 a^{2/3} b}-\frac {\left ((-1)^{2/3} d \sinh \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right )\right ) \int \frac {\cos \left (\frac {\sqrt [6]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-i d x\right )}{\sqrt [3]{a}-\sqrt [3]{-1} \sqrt [3]{b} x} \, dx}{9 \sqrt [3]{a} b^{4/3}} \\ & = -\frac {x \cosh (c+d x)}{3 b \left (a+b x^3\right )}-\frac {\sqrt [3]{-1} \cosh \left (c+\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Chi}\left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right )}{9 a^{2/3} b^{4/3}}+\frac {(-1)^{2/3} \cosh \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Chi}\left (-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right )}{9 a^{2/3} b^{4/3}}+\frac {\cosh \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Chi}\left (\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}+d x\right )}{9 a^{2/3} b^{4/3}}-\frac {d \text {Chi}\left (\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}+d x\right ) \sinh \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{9 \sqrt [3]{a} b^{5/3}}-\frac {(-1)^{2/3} d \text {Chi}\left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right ) \sinh \left (c+\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{9 \sqrt [3]{a} b^{5/3}}+\frac {\sqrt [3]{-1} d \text {Chi}\left (-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right ) \sinh \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{9 \sqrt [3]{a} b^{5/3}}+\frac {(-1)^{2/3} d \cosh \left (c+\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Shi}\left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right )}{9 \sqrt [3]{a} b^{5/3}}+\frac {\sqrt [3]{-1} \sinh \left (c+\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Shi}\left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right )}{9 a^{2/3} b^{4/3}}-\frac {d \cosh \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Shi}\left (\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}+d x\right )}{9 \sqrt [3]{a} b^{5/3}}+\frac {\sinh \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Shi}\left (\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}+d x\right )}{9 a^{2/3} b^{4/3}}+\frac {\sqrt [3]{-1} d \cosh \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Shi}\left (\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}+d x\right )}{9 \sqrt [3]{a} b^{5/3}}+\frac {(-1)^{2/3} \sinh \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Shi}\left (\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}+d x\right )}{9 a^{2/3} b^{4/3}} \\ \end{align*}

Mathematica [C] (verified)

Result contains higher order function than in optimal. Order 9 vs. order 4 in optimal.

Time = 0.18 (sec) , antiderivative size = 363, normalized size of antiderivative = 0.51 \[ \int \frac {x^3 \cosh (c+d x)}{\left (a+b x^3\right )^2} \, dx=\frac {-\frac {6 b x \cosh (c+d x)}{a+b x^3}-\text {RootSum}\left [a+b \text {$\#$1}^3\&,\frac {-\cosh (c+d \text {$\#$1}) \text {Chi}(d (x-\text {$\#$1}))+\text {Chi}(d (x-\text {$\#$1})) \sinh (c+d \text {$\#$1})+\cosh (c+d \text {$\#$1}) \text {Shi}(d (x-\text {$\#$1}))-\sinh (c+d \text {$\#$1}) \text {Shi}(d (x-\text {$\#$1}))+d \cosh (c+d \text {$\#$1}) \text {Chi}(d (x-\text {$\#$1})) \text {$\#$1}-d \text {Chi}(d (x-\text {$\#$1})) \sinh (c+d \text {$\#$1}) \text {$\#$1}-d \cosh (c+d \text {$\#$1}) \text {Shi}(d (x-\text {$\#$1})) \text {$\#$1}+d \sinh (c+d \text {$\#$1}) \text {Shi}(d (x-\text {$\#$1})) \text {$\#$1}}{\text {$\#$1}^2}\&\right ]+\text {RootSum}\left [a+b \text {$\#$1}^3\&,\frac {\cosh (c+d \text {$\#$1}) \text {Chi}(d (x-\text {$\#$1}))+\text {Chi}(d (x-\text {$\#$1})) \sinh (c+d \text {$\#$1})+\cosh (c+d \text {$\#$1}) \text {Shi}(d (x-\text {$\#$1}))+\sinh (c+d \text {$\#$1}) \text {Shi}(d (x-\text {$\#$1}))+d \cosh (c+d \text {$\#$1}) \text {Chi}(d (x-\text {$\#$1})) \text {$\#$1}+d \text {Chi}(d (x-\text {$\#$1})) \sinh (c+d \text {$\#$1}) \text {$\#$1}+d \cosh (c+d \text {$\#$1}) \text {Shi}(d (x-\text {$\#$1})) \text {$\#$1}+d \sinh (c+d \text {$\#$1}) \text {Shi}(d (x-\text {$\#$1})) \text {$\#$1}}{\text {$\#$1}^2}\&\right ]}{18 b^2} \]

[In]

Integrate[(x^3*Cosh[c + d*x])/(a + b*x^3)^2,x]

[Out]

((-6*b*x*Cosh[c + d*x])/(a + b*x^3) - RootSum[a + b*#1^3 & , (-(Cosh[c + d*#1]*CoshIntegral[d*(x - #1)]) + Cos
hIntegral[d*(x - #1)]*Sinh[c + d*#1] + Cosh[c + d*#1]*SinhIntegral[d*(x - #1)] - Sinh[c + d*#1]*SinhIntegral[d
*(x - #1)] + d*Cosh[c + d*#1]*CoshIntegral[d*(x - #1)]*#1 - d*CoshIntegral[d*(x - #1)]*Sinh[c + d*#1]*#1 - d*C
osh[c + d*#1]*SinhIntegral[d*(x - #1)]*#1 + d*Sinh[c + d*#1]*SinhIntegral[d*(x - #1)]*#1)/#1^2 & ] + RootSum[a
 + b*#1^3 & , (Cosh[c + d*#1]*CoshIntegral[d*(x - #1)] + CoshIntegral[d*(x - #1)]*Sinh[c + d*#1] + Cosh[c + d*
#1]*SinhIntegral[d*(x - #1)] + Sinh[c + d*#1]*SinhIntegral[d*(x - #1)] + d*Cosh[c + d*#1]*CoshIntegral[d*(x -
#1)]*#1 + d*CoshIntegral[d*(x - #1)]*Sinh[c + d*#1]*#1 + d*Cosh[c + d*#1]*SinhIntegral[d*(x - #1)]*#1 + d*Sinh
[c + d*#1]*SinhIntegral[d*(x - #1)]*#1)/#1^2 & ])/(18*b^2)

Maple [C] (warning: unable to verify)

Result contains higher order function than in optimal. Order 9 vs. order 4.

Time = 0.42 (sec) , antiderivative size = 877, normalized size of antiderivative = 1.22

method result size
risch \(\text {Expression too large to display}\) \(877\)

[In]

int(x^3*cosh(d*x+c)/(b*x^3+a)^2,x,method=_RETURNVERBOSE)

[Out]

-1/6*d^3*exp(-d*x-c)/b/(b*d^3*x^3+a*d^3)*x-1/18/d/a/b^2*sum((3*_R2^2*b*c^2-_R2*a*d^3-5*_R2*b*c^3-2*a*c*d^3+2*b
*c^4+3*_R2*b*c^2+a*d^3-b*c^3)/(_R2^2-2*_R2*c+c^2)*exp(-_R2)*Ei(1,d*x-_R2+c),_R2=RootOf(_Z^3*b-3*_Z^2*b*c+3*_Z*
b*c^2+a*d^3-b*c^3))+1/18/d*c^3/a/b*sum((_R2-c+2)/(_R2^2-2*_R2*c+c^2)*exp(-_R2)*Ei(1,d*x-_R2+c),_R2=RootOf(_Z^3
*b-3*_Z^2*b*c+3*_Z*b*c^2+a*d^3-b*c^3))-1/6/d*c^2/a/b*sum((_R2^2-_R2*c+_R2+c)/(_R2^2-2*_R2*c+c^2)*exp(-_R2)*Ei(
1,d*x-_R2+c),_R2=RootOf(_Z^3*b-3*_Z^2*b*c+3*_Z*b*c^2+a*d^3-b*c^3))+1/6/d*c/a/b^2*sum((2*_R2^2*b*c-3*_R2*b*c^2-
a*d^3+b*c^3+2*_R2*b*c)/(_R2^2-2*_R2*c+c^2)*exp(-_R2)*Ei(1,d*x-_R2+c),_R2=RootOf(_Z^3*b-3*_Z^2*b*c+3*_Z*b*c^2+a
*d^3-b*c^3))-1/6*d^3*exp(d*x+c)/b/(b*d^3*x^3+a*d^3)*x+1/18/d/a/b^2*sum((3*_R2^2*b*c^2-_R2*a*d^3-5*_R2*b*c^3-2*
a*c*d^3+2*b*c^4-3*_R2*b*c^2-a*d^3+b*c^3)/(_R2^2-2*_R2*c+c^2)*exp(_R2)*Ei(1,-d*x+_R2-c),_R2=RootOf(_Z^3*b-3*_Z^
2*b*c+3*_Z*b*c^2+a*d^3-b*c^3))-1/18/d*c^3/a/b*sum((_R2-c-2)/(_R2^2-2*_R2*c+c^2)*exp(_R2)*Ei(1,-d*x+_R2-c),_R2=
RootOf(_Z^3*b-3*_Z^2*b*c+3*_Z*b*c^2+a*d^3-b*c^3))+1/6/d*c^2/a/b*sum((_R2^2-_R2*c-_R2-c)/(_R2^2-2*_R2*c+c^2)*ex
p(_R2)*Ei(1,-d*x+_R2-c),_R2=RootOf(_Z^3*b-3*_Z^2*b*c+3*_Z*b*c^2+a*d^3-b*c^3))-1/6/d*c/a/b^2*sum((2*_R2^2*b*c-3
*_R2*b*c^2-a*d^3+b*c^3-2*_R2*b*c)/(_R2^2-2*_R2*c+c^2)*exp(_R2)*Ei(1,-d*x+_R2-c),_R2=RootOf(_Z^3*b-3*_Z^2*b*c+3
*_Z*b*c^2+a*d^3-b*c^3))

Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 2050 vs. \(2 (500) = 1000\).

Time = 0.30 (sec) , antiderivative size = 2050, normalized size of antiderivative = 2.86 \[ \int \frac {x^3 \cosh (c+d x)}{\left (a+b x^3\right )^2} \, dx=\text {Too large to display} \]

[In]

integrate(x^3*cosh(d*x+c)/(b*x^3+a)^2,x, algorithm="fricas")

[Out]

-1/36*(12*a*d*x*cosh(d*x + c) - ((a*d^3/b)^(2/3)*((b*x^3 - sqrt(-3)*(b*x^3 + a) + a)*cosh(d*x + c)^2 - (b*x^3
- sqrt(-3)*(b*x^3 + a) + a)*sinh(d*x + c)^2) - (a*d^3/b)^(1/3)*((b*x^3 + sqrt(-3)*(b*x^3 + a) + a)*cosh(d*x +
c)^2 - (b*x^3 + sqrt(-3)*(b*x^3 + a) + a)*sinh(d*x + c)^2))*Ei(d*x - 1/2*(a*d^3/b)^(1/3)*(sqrt(-3) + 1))*cosh(
1/2*(a*d^3/b)^(1/3)*(sqrt(-3) + 1) + c) + ((-a*d^3/b)^(2/3)*((b*x^3 - sqrt(-3)*(b*x^3 + a) + a)*cosh(d*x + c)^
2 - (b*x^3 - sqrt(-3)*(b*x^3 + a) + a)*sinh(d*x + c)^2) - (-a*d^3/b)^(1/3)*((b*x^3 + sqrt(-3)*(b*x^3 + a) + a)
*cosh(d*x + c)^2 - (b*x^3 + sqrt(-3)*(b*x^3 + a) + a)*sinh(d*x + c)^2))*Ei(-d*x - 1/2*(-a*d^3/b)^(1/3)*(sqrt(-
3) + 1))*cosh(1/2*(-a*d^3/b)^(1/3)*(sqrt(-3) + 1) - c) - ((a*d^3/b)^(2/3)*((b*x^3 + sqrt(-3)*(b*x^3 + a) + a)*
cosh(d*x + c)^2 - (b*x^3 + sqrt(-3)*(b*x^3 + a) + a)*sinh(d*x + c)^2) - (a*d^3/b)^(1/3)*((b*x^3 - sqrt(-3)*(b*
x^3 + a) + a)*cosh(d*x + c)^2 - (b*x^3 - sqrt(-3)*(b*x^3 + a) + a)*sinh(d*x + c)^2))*Ei(d*x + 1/2*(a*d^3/b)^(1
/3)*(sqrt(-3) - 1))*cosh(1/2*(a*d^3/b)^(1/3)*(sqrt(-3) - 1) - c) + ((-a*d^3/b)^(2/3)*((b*x^3 + sqrt(-3)*(b*x^3
 + a) + a)*cosh(d*x + c)^2 - (b*x^3 + sqrt(-3)*(b*x^3 + a) + a)*sinh(d*x + c)^2) - (-a*d^3/b)^(1/3)*((b*x^3 -
sqrt(-3)*(b*x^3 + a) + a)*cosh(d*x + c)^2 - (b*x^3 - sqrt(-3)*(b*x^3 + a) + a)*sinh(d*x + c)^2))*Ei(-d*x + 1/2
*(-a*d^3/b)^(1/3)*(sqrt(-3) - 1))*cosh(1/2*(-a*d^3/b)^(1/3)*(sqrt(-3) - 1) + c) - 2*((-a*d^3/b)^(2/3)*((b*x^3
+ a)*cosh(d*x + c)^2 - (b*x^3 + a)*sinh(d*x + c)^2) - (-a*d^3/b)^(1/3)*((b*x^3 + a)*cosh(d*x + c)^2 - (b*x^3 +
 a)*sinh(d*x + c)^2))*Ei(-d*x + (-a*d^3/b)^(1/3))*cosh(c + (-a*d^3/b)^(1/3)) + 2*((a*d^3/b)^(2/3)*((b*x^3 + a)
*cosh(d*x + c)^2 - (b*x^3 + a)*sinh(d*x + c)^2) - (a*d^3/b)^(1/3)*((b*x^3 + a)*cosh(d*x + c)^2 - (b*x^3 + a)*s
inh(d*x + c)^2))*Ei(d*x + (a*d^3/b)^(1/3))*cosh(-c + (a*d^3/b)^(1/3)) - ((a*d^3/b)^(2/3)*((b*x^3 - sqrt(-3)*(b
*x^3 + a) + a)*cosh(d*x + c)^2 - (b*x^3 - sqrt(-3)*(b*x^3 + a) + a)*sinh(d*x + c)^2) - (a*d^3/b)^(1/3)*((b*x^3
 + sqrt(-3)*(b*x^3 + a) + a)*cosh(d*x + c)^2 - (b*x^3 + sqrt(-3)*(b*x^3 + a) + a)*sinh(d*x + c)^2))*Ei(d*x - 1
/2*(a*d^3/b)^(1/3)*(sqrt(-3) + 1))*sinh(1/2*(a*d^3/b)^(1/3)*(sqrt(-3) + 1) + c) + ((-a*d^3/b)^(2/3)*((b*x^3 -
sqrt(-3)*(b*x^3 + a) + a)*cosh(d*x + c)^2 - (b*x^3 - sqrt(-3)*(b*x^3 + a) + a)*sinh(d*x + c)^2) - (-a*d^3/b)^(
1/3)*((b*x^3 + sqrt(-3)*(b*x^3 + a) + a)*cosh(d*x + c)^2 - (b*x^3 + sqrt(-3)*(b*x^3 + a) + a)*sinh(d*x + c)^2)
)*Ei(-d*x - 1/2*(-a*d^3/b)^(1/3)*(sqrt(-3) + 1))*sinh(1/2*(-a*d^3/b)^(1/3)*(sqrt(-3) + 1) - c) + ((a*d^3/b)^(2
/3)*((b*x^3 + sqrt(-3)*(b*x^3 + a) + a)*cosh(d*x + c)^2 - (b*x^3 + sqrt(-3)*(b*x^3 + a) + a)*sinh(d*x + c)^2)
- (a*d^3/b)^(1/3)*((b*x^3 - sqrt(-3)*(b*x^3 + a) + a)*cosh(d*x + c)^2 - (b*x^3 - sqrt(-3)*(b*x^3 + a) + a)*sin
h(d*x + c)^2))*Ei(d*x + 1/2*(a*d^3/b)^(1/3)*(sqrt(-3) - 1))*sinh(1/2*(a*d^3/b)^(1/3)*(sqrt(-3) - 1) - c) - ((-
a*d^3/b)^(2/3)*((b*x^3 + sqrt(-3)*(b*x^3 + a) + a)*cosh(d*x + c)^2 - (b*x^3 + sqrt(-3)*(b*x^3 + a) + a)*sinh(d
*x + c)^2) - (-a*d^3/b)^(1/3)*((b*x^3 - sqrt(-3)*(b*x^3 + a) + a)*cosh(d*x + c)^2 - (b*x^3 - sqrt(-3)*(b*x^3 +
 a) + a)*sinh(d*x + c)^2))*Ei(-d*x + 1/2*(-a*d^3/b)^(1/3)*(sqrt(-3) - 1))*sinh(1/2*(-a*d^3/b)^(1/3)*(sqrt(-3)
- 1) + c) + 2*((-a*d^3/b)^(2/3)*((b*x^3 + a)*cosh(d*x + c)^2 - (b*x^3 + a)*sinh(d*x + c)^2) - (-a*d^3/b)^(1/3)
*((b*x^3 + a)*cosh(d*x + c)^2 - (b*x^3 + a)*sinh(d*x + c)^2))*Ei(-d*x + (-a*d^3/b)^(1/3))*sinh(c + (-a*d^3/b)^
(1/3)) - 2*((a*d^3/b)^(2/3)*((b*x^3 + a)*cosh(d*x + c)^2 - (b*x^3 + a)*sinh(d*x + c)^2) - (a*d^3/b)^(1/3)*((b*
x^3 + a)*cosh(d*x + c)^2 - (b*x^3 + a)*sinh(d*x + c)^2))*Ei(d*x + (a*d^3/b)^(1/3))*sinh(-c + (a*d^3/b)^(1/3)))
/((a*b^2*d*x^3 + a^2*b*d)*cosh(d*x + c)^2 - (a*b^2*d*x^3 + a^2*b*d)*sinh(d*x + c)^2)

Sympy [F(-1)]

Timed out. \[ \int \frac {x^3 \cosh (c+d x)}{\left (a+b x^3\right )^2} \, dx=\text {Timed out} \]

[In]

integrate(x**3*cosh(d*x+c)/(b*x**3+a)**2,x)

[Out]

Timed out

Maxima [F]

\[ \int \frac {x^3 \cosh (c+d x)}{\left (a+b x^3\right )^2} \, dx=\int { \frac {x^{3} \cosh \left (d x + c\right )}{{\left (b x^{3} + a\right )}^{2}} \,d x } \]

[In]

integrate(x^3*cosh(d*x+c)/(b*x^3+a)^2,x, algorithm="maxima")

[Out]

1/2*((d^2*x^3*e^(2*c) + 3*d*x^2*e^(2*c) + 12*x*e^(2*c))*e^(d*x) - (d^2*x^3 - 3*d*x^2 + 12*x)*e^(-d*x))/(b^2*d^
3*x^6*e^c + 2*a*b*d^3*x^3*e^c + a^2*d^3*e^c) + 1/2*integrate(-6*(a*d^2*x^2*e^c - 10*b*x^3*e^c + 3*a*d*x*e^c +
2*a*e^c)*e^(d*x)/(b^3*d^3*x^9 + 3*a*b^2*d^3*x^6 + 3*a^2*b*d^3*x^3 + a^3*d^3), x) - 1/2*integrate(-6*(a*d^2*x^2
 - 10*b*x^3 - 3*a*d*x + 2*a)*e^(-d*x)/(b^3*d^3*x^9*e^c + 3*a*b^2*d^3*x^6*e^c + 3*a^2*b*d^3*x^3*e^c + a^3*d^3*e
^c), x)

Giac [F]

\[ \int \frac {x^3 \cosh (c+d x)}{\left (a+b x^3\right )^2} \, dx=\int { \frac {x^{3} \cosh \left (d x + c\right )}{{\left (b x^{3} + a\right )}^{2}} \,d x } \]

[In]

integrate(x^3*cosh(d*x+c)/(b*x^3+a)^2,x, algorithm="giac")

[Out]

integrate(x^3*cosh(d*x + c)/(b*x^3 + a)^2, x)

Mupad [F(-1)]

Timed out. \[ \int \frac {x^3 \cosh (c+d x)}{\left (a+b x^3\right )^2} \, dx=\int \frac {x^3\,\mathrm {cosh}\left (c+d\,x\right )}{{\left (b\,x^3+a\right )}^2} \,d x \]

[In]

int((x^3*cosh(c + d*x))/(a + b*x^3)^2,x)

[Out]

int((x^3*cosh(c + d*x))/(a + b*x^3)^2, x)